Techniques for howling detection

ABSTRACT

A howling detector is described that is configured to receive an input signal and to determine measure of the linearity of a logarithmic representation of the energy of the input signal. In some examples, this triggers gain adjustment (e.g. of a noise control unit) and, in some further examples, the amount of the gain adjustment may be based on an estimation of the maximum stable gain of a noise control unit.

The present application is a 371 of International Patent Application No. PCT/GB2019/051467, filed May 29, 2019, which claims priority to U.S. patent application Ser. No. 16/004,900, filed Jun. 11, 2018, now U.S. Pat. No. 10,681,458, issued Jun. 9, 2020, and United Kingdom Patent Application Serial No. 1811398.5, filed Jul. 12, 2018, each of which is incorporated by reference herein in its entirety.

TECHNICAL FIELD

Embodiments described herein relate to techniques for active noise cancellation. In particular, embodiments described herein relate to techniques for suppressing an acoustic feedback signal.

BACKGROUND

Many devices include loudspeakers, which are used to play sounds to a user of the device, based on an input signal. For example, the input signal may be derived from a signal that has been received by the device over a communications link, in the case of a phone call or the like, or may be derived from stored data, in the case of music or speech playback.

As wireless communication devices, Mp3 players and other devices for audio playback move even further into everyday use, features like noise cancellation become more important to help ensure higher-quality audio playback and phone calls.

Noise cancellation embraces a number of different approaches to eliminating unwanted noise in order to enhance the listening experience of a user. Active noise cancellation (ANC) or noise control refers to a method of reducing noise by the addition of anti-noise—i.e. a phase inverted noise signal—which destructively interferes with the noise. This is generally achieved by using a reference microphone to sense environmental or ambient noise and by deriving an anti-noise signal that is emitted by a speaker in order to cancel or at least control the noise. As will be appreciated by those skilled in the art active noise control can be achieved with analogue filters or digital filters, and is generally differentiated by architecture: feed-forward cancellation, feedback cancellation or hybrid cancellation.

FIG. 1 provides a simplified illustration of a feedforward ANC system. As illustrated in FIG. 1 , a reference microphone 10 detects incident ambient sounds—or noise—and generates an input signal x(n) for an ANC circuit 20. The ANC circuit 20 processes the signal in order to derive a control signal y(n) which is passed to the loudspeaker transducer 30 and is emitted by the loudspeaker 30 as anti-noise. Thus, the ANC circuit may be considered to comprise a filter having a transfer function H_(nc) which inversely models the noise signal for generating the required control sound signal. An error microphone (not shown) is typically provided to measure the error between the noise signal and the anti-noise signal in order that the transfer function H_(nc) or the respective gain of the ANC circuit may be adapted.

As illustrated in FIG. 2 , it will be appreciated that the anti-noise signal will not only propagate on a path towards a user's ear H_(de) (where d denotes the driver and e denotes the ear), but may also propagate on a leakage path, or feedback path H_(dm) (where d denotes the driver/loudspeaker and m denotes the microphone), towards the reference microphone. This is known as acoustic feedback and results in a corrupted reference signal u(n). Thus, the reference signal will additionally contain the acoustic feedback signal that is sensed by the reference microphone. When an acoustic control system has a feedback path the leakage often causes unstable behaviour called howling which results in an audible feedback tone.

It will therefore be appreciated that the stability of a noise control system will be significantly influenced by the feedback signal and will depend on the transfer characteristics of an acoustic feedback path H_(dm) between the speaker and the reference microphone. A similar problem can arise when the speech captured by the voice microphone leaks to the speaker (driver). The problem of acoustic feedback is particularly an issue in the case of a mobile communication device, such as a mobile phone, due to the close proximity between the reference microphone and the speaker.

The frequency of the feedback tone depends on Hdm in conjunction with Hnc. Since both Hdm and Hnc can change, the tone frequency can change respectively across a wide range of frequencies. The level of the feedback tone may rise quickly as the energy of the reference signal rises exponentially. Thus, the feedback tone is unpleasant and potentially damaging to the ear. There is therefore a need to try to manage and/or suppress the occurrence of howling.

A number of techniques have been proposed that seek to detect the occurrence of howling in order to enable the gain of the circuit to be automatically adjusted. However, the previously considered techniques suffer from a number of disadvantages, including high latency, the need for relatively complex frequency domain processing and problems associated with the occurrence of false positive detection of howling. For example, a previously proposed method involves a maximum peak detection method which involves performing a linear search for a sustained peak in the energy of the signal across multiple consecutive frames. As a consequence of the need to locate and track several frames of the signal, the system introduces a latency between the initial occurrence of howling and the detection of howling. Consequently, any subsequent measures taken to mitigate the feedback tone take place after a certain delay, and potentially after the sound has risen to an audible level.

The present examples are concerned with techniques for detecting howling, in particular to howling detection techniques which alleviate one or more of the problems associated with previously proposed howling detection methods.

According to one example of a first aspect of the present invention there is provided a howling detector (or howling detection unit) configured to receive an input signal and to determine a measure of the linearity of a logarithmic representation of the energy of the input signal. As will be explained below, according to one or more examples herein the measure of linearity (or linearity metric) of the logarithm of the energy may constitute an indication that howling is likely, imminent, or occurring and therefore the howling detector of some examples herein is configured to detect the presence of howling. In some examples, the input signal may comprise a feedback signal.

According to one or more examples, the howling detector may be configured to issue a command to a gain adjuster if the measure of linearity exceeds a predetermined threshold. The command issued to the gain adjuster may comprise instructions that, when executed by the gain adjuster, cause the gain adjuster to reduce the gain. The command issued to the gain adjuster may comprise instructions that, when executed by the gain adjuster, cause the gain adjuster to adjust the gain, for example by a fixed amount. In this way, once the howling detector detects howling (or imminent or likely howling), it may issue a command to a gain adjuster to suppress the howling. As the howling may be configured to automatically issue the command to the gain adjuster the howling detector is therefore configured to actively suppress, or cancel, the effects of howling. It will be appreciated that, in some examples, the measure of linearity exceeding a predetermined threshold may be determined by a parameter, or phasemeters, falling below a predetermined threshold. In other words, a parameter (examples will be provided below and include a residual value or a derivative) falling below a threshold may indicate that the measure of linearity of the logarithm of energy exceed a threshold. In other words, certain parameters may decrease the more linear the logarithm of the energy becomes.

According to one or more examples, the howling detector may be configured to issue a command to a gain adjuster if the measure of linearity falls below a predetermined threshold, depending on the measure of linearity as will be explained below.

Determining a measure of the linearity of the logarithmic representation of the energy of the input signal may, in some examples, comprise computing a difference between the logarithm of the energy of the signal and a trend line, the trend line being a straight line representation of the logarithm of the energy of the signal. In this way, the measure of linearity may be the difference between the logarithm of the signal energy and a straight line representation of the logarithm of the signal energy, and therefore when the logarithm of the energy tends towards linearity this difference will become smaller. This difference may be measured by the residuals of the trend line.

Accordingly, in one example, the howling detector may be configured to issue command to a gain adjuster if the average of the residuals of the trend line are below a threshold. In this way, when the residuals of the trend line are below a threshold this may be an indication that the logarithm of the energy tends sufficiently towards linearity so as to infer not just the presence of howling, but when howling is likely or imminent (depending on the threshold value). Issuing the command to the gain adjuster when the residuals fall below a threshold means that the howling detector may instruct the gain adjuster to reduce the gain when howling is likely or imminent, and therefore, in some examples, before howling has actually occurred.

According to one or more examples, coefficients of the trend line are derived by fitting a line to the logarithmic representation of the energy of the input signal. Fitting a line to a logarithmic representation of the energy of the input signal may comprise performing a least squares computation which minimises a difference between the trend line and the logarithmic representation of the energy of the input signal.

The measure of linearity may be determined according to a goodness of fit measure of the computed trend line to the logarithm of the energy of the input signal over P samples. In one example, the goodness of fit measure r(n) may be represented by:

${{r(n)} = \sqrt{\sum\limits_{i = 0}^{P - 1}\left( {{y\left( {n - i} \right)} - \left( {{{\hat{\beta}}_{0}(n)} + {\left( {P - i} \right) \cdot {{\hat{\beta}}_{1}(n)}}} \right)} \right)^{2}}},$ where y(n) is the log energy of the input signal, {circumflex over (β)}₀(n) is the bias of the trend line and {circumflex over (β)}₁(n) is the slope of the trend line.

The howling detector may be configured to issue the command to the gain adjuster when it is detected that the slope {circumflex over (β)}₁(n) is converging to the value

$\frac{1}{P - 1}.$ In this example, convergence of the slope may sufficiently indicate that the logarithm of the signal energy is linear, or sufficiently tending towards linearity. Convergence of the slope may be measured by the slope being within a confidence interval and, accordingly, detecting that the slope is converging may comprise determining whether the slope lies within a confidence interval, for example the interval:

${{\left( \frac{1}{P - 1} \right) \cdot \left( {1 - \alpha} \right)} \leq {{{\overset{\_}{\beta}}_{1}\left( \frac{1}{P - 1} \right)} \cdot \left( {1 + \alpha} \right)}},$ where α is a parameter.

The howling detector may be configured to determine a derivative of the logarithm of the energy. The howling detector may be configured to issue the command to the gain adjuster when the derivative of the logarithm of the energy tends towards zero or is less than a threshold.

According to an example of a second aspect there is provided circuitry for receiving an input signal and determining a measure of the linearity of the logarithmic representation of the energy of the input signal. The input signal may comprise a feedback signal.

The circuitry may, in some examples, comprise a gain adjuster configured to reduce the gain if the measure of linearity exceeds a certain threshold. In these examples, the circuitry may be configured to control the gain adjuster to reduce the gain if the measure of linearity exceeds a predetermined threshold. The gain adjuster may be configured to adjust the gain by an amount that relates to a difference between a measure of the linearity of the plot of the energy of the input signal in the log domain and a target linearity metric. For example, the gain adjuster may be configured to adjust the gain by a fixed amount, or by an amount which is proportional to the slope of the trend line (as described above).

The circuitry may comprise a filterbank configured to split the input signal into a plurality of frequency bands, wherein the circuitry is operable to determine a measure of the linearity of the logarithmic representation of the energy for each frequency band.

The circuitry may comprise the howling detector as described above. In this example, the howling detector may be configured to issue a command to the gain adjuster if the measure of linearity exceeds a predetermined threshold.

The circuitry may comprise a noise control module, or a noise control system (e.g. a noise reduction system) such as active noise control module or active noise control system. The circuitry may comprise processing circuitry for one of the aforementioned modules or systems.

According to one example of a third aspect there is provided a method comprising: receiving an input signal, determining a measure of the linearity of the logarithmic representation of the energy of the input signal.

The method may further comprise issuing a command to a gain adjuster if the measure of the linearity exceeds a predetermined threshold.

Determining a measure of the linearity of the logarithmic representation of the energy of the input signal may comprise computing a difference between the logarithm of the energy of the signal and a trend line, the trend line being a straight line representation of the logarithm of the energy of the signal, and the method may further comprise issuing the command to the gain adjuster if the average of the residuals of the trend line are below a threshold.

Determining a measure of the linearity of the logarithmic representation of the energy of the input signal may comprise calculating a goodness of fit measure of a computed trend line being a straight line representation of the logarithm of the energy of the signal to the logarithm of the energy of the input signal over P samples, wherein the goodness of fit measure r(n) is represented by:

${{r(n)} = \sqrt{\sum\limits_{i = 0}^{P - 1}\left( {{y\left( {n - i} \right)} - \left( {{{\hat{\beta}}_{0}(n)} + {\left( {P - i} \right) \cdot {{\hat{\beta}}_{1}(n)}}} \right)} \right)^{2}}},$ where y(n) is the log energy of the input signal, {circumflex over (β)}₀(n) is the bias of the trend line and {circumflex over (β)}₁(n) is the slope of the trend line, and the method may further comprise issuing the command when it is detected that the slope {circumflex over (β)}₁(n) is converging to the value

$\frac{1}{P - 1}.$

Determining a measure of the linearity of the logarithmic representation of the energy of the input signal may comprise determining a derivative of the logarithm of the energy, and the method may further comprise issuing the command to the gain adjuster when the derivative of the logarithm of the energy tends towards zero or is less than a threshold.

According to an example of a fourth aspect there is provided a processing module for a noise control circuit, the processing module comprising:

-   -   a howling detector configured to receive an input signal and to         determine a linearity metric based on the input signal, the         linearity metric comprising a measure of the linearity of a         logarithmic representation of the energy of the input signal;         and     -   a gain adjuster configured to adjust the gain of a noise control         unit.

The linearity metric—or log-linearity metric—may be derived by computing a difference between the logarithm of the energy of the signal and a trend line, the trend line being a straight line representation of the logarithm of the energy of the signal. Coefficients of the trend line may be derived by a fitting a line to the logarithmic representation of the energy of the input signal. The process of fitting a line to the logarithmic representation of the energy of the input signal may comprise performing a least squares computation which minimises a difference between the trend line and the logarithmic representation of the energy of the input signal.

According to one or more examples the linearity metric is determined according to a goodness of fit measure of the computed trend line to the logarithm of the energy of the input signal over P samples. The goodness of fit measure r(n) may be represented by:

${r(n)} = \sqrt{\sum\limits_{i = 0}^{P - 1}\left( {{y\left( {n - i} \right)} - \left( {{{\hat{\beta}}_{0}(n)} + {\left( {P - i} \right) \cdot {{\hat{\beta}}_{1}(n)}}} \right)} \right)^{2}}$ where y(n) is the log energy of the input signal, {circumflex over (β)}₀(n) is the bias of the trend line and {circumflex over (β)}₁(n) is the slope of the trend line.

According to one or more examples the howling detection unit may be further configured to issue a command to the gain adjuster in order to reduce the gain if the linearity metric exceeds a predetermined threshold. The gain adjuster may be configured to adjust the gain by a fixed amount or by an amount which is proportional to the slope of the trend line.

According to one or more examples the howling detection unit may be further configured to estimate a maximum stable gain of the noise control unit. The maximum stable gain is proportional to the slope of the trend line.

The howling detection unit may be configured to determine, based on the determined linearity metric, if howling is likely or imminent.

According to a further aspect there is provided an audio processing system comprising a processing module according to the fourth aspect. The audio processing system may further comprise a noise control unit for generating a noise control signal based on a reference input signal which represents a sound detected by a reference microphone. The processing module may be connected to a speaker and wherein the speaker generates an anti-noise signal based on the noise control signal in order to cancel or at least reduce the noise detected by a reference microphone.

According to one example the howling detector may be provided in parallel with the noise control unit.

According to one example the audio processing system may further comprise a filter configured to filter out one or more frequencies or frequency bands of the input signal. The audio processing system may comprise a filterbank configured to split the input signal into a plurality of frequency bands, wherein the howling detector is operable to determine a linearity metric for each frequency band.

According to at least one example of a sixth aspect, there is provided a processing circuit for a noise control module comprising: a gain adjustment mechanism configured to adjust the gain of a noise control circuit if a plot of the energy of an input signal in the log domain becomes linear or tends towards linearity.

The gain of the noise control circuit may be adjusted by an amount that relates to an amount which is proportional to the slope of the trend line.

The processing module may be provided in the form of a monolithic integrated circuit.

According to an example of a further aspect there is provided a device comprising a processing module according to the fourth aspect or a processing circuit according to the sixth aspect, wherein the device comprises a mobile telephone, headphone, acoustic noise cancelling headphones, a smart watch, an audio player, a video player, a mobile computing platform, a games device, a remote controller device, a toy, a machine, or a home automation controller, a domestic appliance or other portable device.

According to a seventh aspect there is provided a method of processing an audio signal comprising:

-   -   determining a linearity metric based on an input signal, the         linearity metric comprising a measure of the linearity of a         logarithmic representation of the energy of the input signal;         and     -   adjusting the gain of the noise control unit if the linearity         metric exceeds a predetermined threshold.

The step of determining a linearity metric may comprise computing a difference between the logarithm of the energy of the signal and a trend line, the trend line being a straight line representation of the logarithm of the energy of the signal.

Aspects of the invention therefore relate to the detection of occurring, imminent, or likely to be imminent, howling and its suppression. For example, some examples relate to the monitoring of parameters that, by exceeding or falling below a pre-determined threshold will indicate a likely to be imminent, or imminent, occurrence of howling. When this occurs a signal may be sent to adjust the gain so that howling can be actively suppressed before howling actually occurs. In this way, some aspects of the invention may be used as part of an active noise reduction, or control, system, that will positively act to suppress howling before it has occurred. Some aspects of the invention may therefore be used in conjunction with an active noise reduction, or control system. However, in, further examples, some aspects may be used as a standalone apparatus to monitor and/or detect the presence of howling.

BRIEF DESCRIPTION OF DRAWINGS

For a better understanding of the present invention and to show how the same may be carried into effect, reference will now be made by way of example to the accompanying drawings in which:

FIG. 1 provides a simplified illustration of a feedforward ANC system;

FIG. 2 provides a further illustration of a feedforward ANC system;

FIG. 3 illustrates a wireless communication device implementing a noise control circuit in accordance with a first example;

FIG. 4 is a block diagram of the audio integrated circuit provided in the wireless communication device depicted in FIG. 3 ;

FIG. 5 illustrates a first example of a processing module;

FIG. 6 illustrates a further example of a processing module according to the present aspects implemented within a noise control module 140;

FIGS. 7 a and 7 b show examples of a line fitting method;

FIG. 8 illustrates a further example of a howling detection circuit implemented within an audio processing circuit;

FIG. 9 shows first and second plots to illustrate the increasing slope in the energy or of an input signal;

FIG. 10 shows a plurality of data samples plotted on a graphical representation of the amplitude of an input signal, as well as a fitted line;

FIG. 11 shows the log energy of the signal output for 3 different cases;

FIG. 12 shows several plots to illustrate the signals used by the proposed method to detect howling;

FIG. 13 illustrates a howling detection process according to one example; and

FIG. 14 illustrates a system incorporating multiband howling detection.

DETAILED DESCRIPTION

The description below sets forth examples according to the present disclosure. Further example embodiments and implementations will be apparent to those having ordinary skill in the art. Further, those having ordinary skill in the art will recognize that various equivalent techniques may be applied in lieu of, or in conjunction with, the examples discussed below, and all such equivalents should be deemed as being encompassed by the present disclosure.

The methods described herein can be implemented in a wide range of devices and systems. However, for ease of explanation of one example, an illustrative example will be described, in which the implementation occurs in a mobile communication device such as a mobile phone.

FIG. 3 illustrates a wireless communication device 100 implementing a noise control unit in accordance with any of the present aspects. The wireless communication device 100 is shown in proximity with a user's ear 50.

The wireless communication device comprises a transducer, such as a speaker 130, which is configured to reproduce distance sounds, such as speech, received by the wireless communication device along with other local audio events such as ringtones, stored audio program material, and other audio effects including a noise control signal. A reference microphone 110 is provided for sensing ambient acoustic events. The wireless communication device further comprises a near-speech microphone 150 which is provided in proximity to a user's mouth to sense sounds, such as speech, generated by the user.

A circuit 125 within the wireless communication device comprises an audio CODEC integrated circuit (IC) 180 that receives the signals from the reference microphone, the near-speech microphone 150 and interfaces with the speaker and other integrated circuits such as a radio frequency (RF) integrated circuit 12 having a wireless telephone transceiver.

FIG. 4 is a block diagram of the audio integrated circuit 180 provided in the wireless communication device depicted in FIG. 3 and illustrates selected units of the integrated circuit. Specifically, the integrated circuit receives an input signal from reference microphone 110 and includes an analog-to-digital converter (ADC) 135 a for generating a digital representation of the input signal x(n) which is passed to a noise control unit 140 according to an example of the present aspects, wherein the noise control circuit is configured to generate a control signal u(n).

The audio integrated circuit comprises a further ADC 135 b for generating a digital representation of the signal generated by the near-speech microphone 150. Combiner 136 may combine audio signals including the noise control signal u(n) (which by convention may have the same polarity as the noise in the reference microphone signal and will therefore be subtracted by the combiner 136), a portion of the near-speech microphone signal to allow a user of the wireless communication device to hear his or her own voice, in addition to downlink speech communication which is received from the radio frequency integrated circuit 115. The digital-to-analogue converter 137 receives the output of the combiner 136, amplifies the output of the DAC (not shown) and passes the resultant signal to the speaker 130.

From consideration of FIG. 3 it will be appreciated that a feedback path H_(dm) arises between the speaker 130 and the reference microphone 110.

FIG. 5 shows a processing module (howling detection system) 200 comprising a howling detector (HD) 150 and a gain adjuster 160 (e.g. amplifier) configured to control the gain of a noise control system comprising a noise control unit 120. The howling detector is configured to determine a linearity metric, or log-linearity metric, which is indicative of a measure of the linearity of the variation (over P samples) in the energy of an input signal u(n) in the logarithm domain. Thus, one or more of the present examples can be considered to rely upon a consideration of the change in the energy of the signal u(n) in the log domain.

According to one or more examples the linearity metric may be derived by computing a difference between the logarithm of the energy of the signal u(n) and a line —referred to herein as a trend line—which is a straight line representation of the logarithm of the energy of the signal u(n). It will be appreciated that the trend line may be derived by fitting a line of best fit to the log-energy data derived from the observed signal samples.

According to one or more examples the linearity metric may be determined according to a goodness of fit measure of the computed trend line to the logarithm of the energy of the input signal u(n) over P samples. The goodness of fit measure may be computed, for example, as the least square L²-norm of the residuals r(n) of the trend line fitted to the observed energy data points. Thus, r(n) may be represented by:

$\begin{matrix} {{r(n)} = \sqrt{\sum\limits_{i = 0}^{P - 1}\left( {{y\left( {n - i} \right)} - \left( {{{\hat{\beta}}_{0}(n)} + {\left( {P - i} \right) \cdot {{\hat{\beta}}_{1}(n)}}} \right)} \right)^{2}}} & {{Eq}.(1)} \end{matrix}$

Where y(n) is the log energy of the signal u(n), {circumflex over (β)}₀(n) is the bias of the trend line and {circumflex over (β)}₁(n) is the slope of the trend line. Thus, equation 1 can be understood to represent the least squared difference between the logarithm of the energy of the signal y(n) and a straight line representation of the logarithm of the energy of the signal u(n). In other words, the linearity metric that is determined by the howling detector 150 according to the present examples may be considered to be r(n).

FIGS. 7 a and 7 b show examples of a line fitting method. In both plots the solid line represents the output log energy y(n) and the dashed line is the computed trend line, or fitted line, having line equation coefficients {circumflex over (β)}(n). In FIG. 7 a a linearity metric derived on the basis of the computed trend line will be relatively low, indicating high residual values r, because the line doesn't closely represent y(n). However, in FIG. 7 b a linearity metric derived on the basis of the computed trend line will be relatively high, indicating low residual values r, because the signal y(n) is closely represented by the fitted line. Thus, it will be appreciated that the linearity metric or goodness of fit measure can be used as an indicator of howling by detecting when the linearity metric exceeds a threshold value (in other words when the residual values r(n) being below a given threshold). It will also be appreciated that as the logarithm of the signal energy tends towards linearity its derivative will tend towards zero, the gradient or slope of straight lines being zero. Therefore, the derivative of the logarithm of the energy may be used as an indicator of howling by detecting when the derivative is below a given threshold (indicating that the line that is the logarithm of the energy is tending towards a zero gradient).

According to one or more examples, if the linearity metric determined by the howling detector 150 exceeds a predetermined threshold this is an indication that howling is occurring.

According to one or more examples the howling detection unit 150 is configured to issue a command to the gain adjuster 160 to reduce the gain of the noise control unit.

According to one or more examples the gain is reduced by a fixed amount in order to effectively reduce or turn off the noise control functionality. Thus, the howling detection circuit may be configured to cause the gain adjuster to drop the gain by a fixed amount that is pre-selected in accordance with the operating conditions of the overall system. For example, the howling detection circuit may be configured to cause the gain adjuster to drop the gain by between 25 and 40 dB, preferably by 35 dB in order to immediately kill the tone, interrupt the function of the noise control unit and alleviate any other audio distortion.

FIG. 6 illustrates a possible implementation of a howling detection circuit, comprising a howling detector 150 and a gain adjuster 160 within a noise control module 140. The noise control module or noise control circuit may be implemented within a monolithic integrated circuit such as the CODEC 180 illustrated in FIG. 4 , or may be implemented within any other audio processing circuit which incorporates a noise control unit for controlling or cancelling noise. It will also be appreciated that the noise control module may be provided in a variety of different devices, not just within a mobile phone.

The noise control module 140 comprises a noise control unit 120 which receives, e.g. via an analogue-to-digital converter (not shown), an input signal x(n) derived from the output signal of a reference microphone 110. The input signal represents ambient noise and other sounds including a feedback signal h_(dm) detected by the reference microphone. The noise control unit 120 is operable to generate a control signal u(n) based on the input signal x(n). The control signal comprises an anti-noise signal—in other words a signal representing the noise but with inverted phase (antiphase)—that is emitted by a speaker 130 of the device in which the noise control circuit is implemented, in order to cancel or at least control or mitigate the level of the noise that is heard by a user.

In this example the control signal forms the input signal u(n) for the howling detection system 200.

FIG. 8 illustrates a further example in which a howling detection circuit 200 may be provided in conjunction with an automatic gain adjustment (AGA) circuit 170 which is operable, based on an input signal derived from an error microphone of the noise control system, to track the energy of the input signal over multiple frames and to make adjustments to the gain of the ANC circuit via the gain adjuster 160. Thus, the gain G may be primarily controlled by the AGA system in order to minimise the residual noise at the speaker 130. The howling detection circuit may be configured to take control of the gain adjuster 160 if howling is detected. If howling occurs, the howling detection circuit is configured to cause the gain adjuster to drop the gain by around 25 to 40 dB, in order to interrupt the operation of the noise control unit 120 and bring a stop to the audio distortion.

One or more examples described herein may advantageously allow howling detection based on a single frame (e.g. 5 msec). Thus, the present examples potentially benefit from a faster detection of howling than previously proposed techniques which require computations to take place over multiple frames. Furthermore, examples of the present aspects are advantageously performed purely in the time domain, circumventing the need to transform the signal to frequency domain, thereby lowering the complexity of the system and easing the memory capacity requirements in comparison to techniques which require frequency domain processing.

A theoretical background to the present examples will now be discussed with reference to FIG. 9 .

When the howling occurs, an exponentially rising sinusoid appears in the signal of interest u(n). The increasing slope in the energy or amplitude of u(n), which is illustrated graphically by the first plot shown in FIG. 9 , is clearly exponentially dependent on the gain. Thus, it will be appreciated that if the logarithm of the energy of the signal is computed, the dependence will become linear, as illustrated in the second plot shown in FIG. 9 . In other words, when howling is present the log energy increases linearly.

Examples described herein are based on monitoring for an exponential increase in the amplitude of the signal, which may be readily recognised from a linear increase in the log-domain (dB), in order to try to detect or predict the occurrence of howling.

The output log energy y(n) may be computed using a smoothing of the instant energy of the input signal u(n). Thus: y(n)=log(a(n)), where a(n)=λ_(y) ·a(n−1)+(1−λy)·u ²(n)  Eq. (2)

The energy a(n) is therefore proportional to the square of the sample.

If we define a straight line where {circumflex over (β)}₀(n) is the bias of the line and {circumflex over (β)}₁(n) is the slope of the line, the line equation term in equation 1 can be written in the vector format as: {circumflex over (β)}₀+(P−i)·{circumflex over (β)}₁=[{circumflex over (β)}₀(n),{circumflex over (β)}₁(n)]^(T)·[1,P−i]={circumflex over (β)}(n)·[1,P−i]  Eq. (3)

Where T indicates matrix transpose operation.

The line equation coefficients {circumflex over (β)}(n) can be computed by minimising a cost function between observed samples and fitted line equation. This can be done for example by solving a least squares problem which tries to find the line {circumflex over (β)}(n) that minimises the Euclidean distance to the observed data y(n).

$\begin{matrix} {{\hat{\beta}(n)} = {{\arg{\min\limits_{\beta}\left( {r(n)}^{2} \right)}} = {\left( {X^{T}X} \right)^{- 1}X^{T}{y(n)}}}} & {{Eq}.(4)} \end{matrix}$

Where:

${X = \begin{bmatrix} 1 & 1 \\ 1 & 2 \\  \vdots & \vdots \\ 1 & P \end{bmatrix}},{{y(n)} = \begin{bmatrix} \begin{matrix} \begin{matrix} {y\left( {n - P + 1} \right)} \\ {y\left( {n - P + 2} \right)} \end{matrix} \\  \vdots  \end{matrix} \\ {y(n)} \end{bmatrix}}$

In order to make this approach robust to different energy levels of the incoming signal u(n), the signal y(n) may be normalised such that it is between 0 and 1.

$\begin{matrix} {{y^{\prime}(n)} = \frac{{y(n)} - {\min\left( {y(n)} \right)}}{{\max\left( {y(n)} \right)} - {\min\left( {y(n)} \right)}}} & {{Eq}.(5)} \end{matrix}$

Due to this normalisation, the slope {circumflex over (β)}₁(n) must converge to

$\frac{1}{P - 1}.$ in the presence of howling and the line coefficients may be computed using:

${\hat{\beta}(n)} = {{\left( {\begin{bmatrix} 1 & 1 & \ldots & 1 \\ 1 & 2 & \ldots & P \end{bmatrix} \cdot \begin{bmatrix} 1 & 1 \\ 1 & 2 \\  \vdots & \vdots \\ 1 & P \end{bmatrix}} \right)^{- 1} \cdot \begin{bmatrix} 1 & 1 & \ldots & 1 \\ 1 & 2 & \ldots & P \end{bmatrix} \cdot \begin{bmatrix} \begin{matrix} \begin{matrix} {y^{\prime}\left( {n - P + 1} \right)} \\ {y^{\prime}\left( {n - P + 2} \right)} \end{matrix} \\  \vdots  \end{matrix} \\ {y^{\prime}(n)} \end{bmatrix}} = {\overset{\sim}{X} \cdot \begin{bmatrix} \begin{matrix} \begin{matrix} {y^{\prime}\left( {n - P + 1} \right)} \\ {y^{\prime}\left( {n - P + 2} \right)} \end{matrix} \\  \vdots  \end{matrix} \\ {y^{\prime}(n)} \end{bmatrix}}}$

Where x may be precomputed to avoid matrix inversion and to reduce the number of cycles. In addition, in order to reduce MIPS, y′(n) can be subsampled to only use a subset of the P samples to solve the least squares problem presented in Eq. 3. In the example illustrated in FIG. 10 the new subset only includes 16 samples represented by each of the circles and obtained from y(n), the samples being equidistant apart starting from first sample and ending with the last sample.

Smoothing may be applied to the estimated slope coefficient {circumflex over (β)}₁(n) and to the residual r(n) r (n)=λ_(r) ·T(n−1)+(1−λ_(r))·r(n)  Eq. (7) {circumflex over (β)}₁(n)=λ_(β)·β ₁(n−1)+(1−λ_(β))·{circumflex over (β)}₁(n)  Eq. (8)

In addition, according to one or more examples the parameter α is used as a confidence interval on the estimated slope coefficient such that howling is only determined to be detected if the slope is between these confidence intervals.

$\begin{matrix} {{\left( \frac{1}{P - 1} \right) \cdot \left( {1 - \alpha} \right)} \leq {\overset{\_}{\beta}}_{1} \leq {\left( \frac{1}{P - 1} \right) \cdot \left( {1 + \alpha} \right)}} & {{Eq}.(9)} \end{matrix}$

According to one or more examples, the predetermined threshold of the linearity metric is set such that a positive indication of howling is determined before log-linearity is actually reached. The positive indication may therefore beneficially provide an advance indication that howling is likely to occur or is imminent. Such examples may be implemented in conjunction with a process which additionally monitors the change in the linearity metric over time. It will be appreciated that it is possible to detect when the linearity metric increases—i.e. when the logarithm of the energy of the signal u(n) tends towards, or approaches, a certain degree of linearity. Thus, according to one or more examples, once the linearity metric meets a predetermined threshold, a further detection process is carried out to ascertain if the linearity metric is tending towards or is approaching a value which indicates a high degree of linearity—i.e. that the linearity is increasing. Specifically, a reduction in the value of the residuals r(n) over time will indicate that the system is tending towards howling.

According to one or more examples, the predetermined threshold of the measure of linearity may be a maximum threshold. In other words, a positive indication of howling may be determined when the measure of linearity falls below a predetermined threshold. In these examples, the positive indication may still provide an advance indication that howling is likely to occur or is imminent. Such examples may be implemented in conjunction with a process which monitors the rate of change in the logarithm of the signal energy. It will be appreciated that when this derivative is zero then the logarithm is perfectly linear, and so the presence of howling may therefore be inferred by monitoring the derivative and determining when it is tending towards zero. In other words, a process may monitor Dy(n)=D log(u(n)), D being the general term used to denote the derivative, which may for example (with reference to FIGS. 7, 9 and 10 ) may be the derivative with respect to samples, and may provide an advance indication that howling is likely to occur when Dy(n)=0. In one example an advance indication that howling is likely to occur may be provided when the derivative of the logarithm of the energy is below a threshold, in other words when Dy(n)<T, T representing a threshold that is approximately zero. In this way when the derivative of the logarithm of the energy is tending towards zero, as indicated by it being below the threshold T, this may provide an advance indication that howling is likely to occur.

Thus, a further advantage of monitoring the linearity of the log-energy is that it is possible to observe when the system is tending towards howling (as indicated e.g. by a goodness of fit measure). In other words, according to one or more examples, it is possible to observe when the system is starting to become unstable, rather than relying upon the detection of a characteristic of howling after it has started. Thus, according to one or more examples, it is possible to determine when howling is likely, or imminent. This will be prior to the actual occurrence of any audible howling.

The linearity metric determined by the howling detector 150 can therefore be considered to be a representation or estimation of the stability of the noise control system. If the system is tending towards howling, which is an unstable condition, adjustment may be beneficially made to the gain G of the circuit in order to restore the noise control system to a more stable state. According to one or more examples, the gain is reduced by an amount which is related to a measure of the degree of instability. Thus, according to one or more examples the gain may be reduced by an amount which depends on how unstable the system is. Thus, the howling detector may be configured to estimate a required amount of gain reduction that is proportional to the degree of instability. This may, for example, be derived from the slope of the un-normalised log-energy. This slope, when howling is present, is proportional to the gain that the system needs to reduce to keep the system stable (no howling).

Thus, the amount by which the gain needs to be adjusted can be considered to be a maximum stable gain (MSG) of the system. The maximum stable gain (MSG) can be considered to be a measure that shows if the system is stable (MSG negative) or if it is unstable (MSG positive). This is actually the minimum extra gain that the system needs to decrease to become stable.

The information from the slope can be also used to adjust the location of poles and zeros of the feed-forward noise cancelling filter.

An advantage of determining an amount of gain adjustment required based on the slope of the fitted trend line that it becomes possible to keep the maximum gain in the ANC without making the system howl. Thus, some degree of noise control can still take place even though the system is in or tending towards an unstable state. In contrast, if the gain is reduced by a fixed amount (35 dB), this noise control functionality is effectively prevented.

FIG. 11 shows the log energy of the signal output for three different cases (values of MSG) where the system becomes unstable after the sample 48000. From this it is apparent that the slope of the fitted line is directly related to the desired MSG value. Thus, according to one or more examples, it is possible to determine the minimum gain reduction necessary.

FIG. 12 shows several plots relating to a test signal and signals used by the proposed method to detect howling. The first plot shows the test signal in the time domain over 3000 samples (approximately 62 ms). The second plot displays the computed log energy from the signal shown first plot. The third and fourth plots show the two metrics used to decide if howling is present. In this example, the system becomes unstable around sample 2300. At this point, the log energy y(n) starts to clearly grow at a linear rate which makes the residual smaller since y(n) is closer to a line. Furthermore, the slope of this log energy y(n) (last figure) is also converging to the expected value plotted as a red line in the figure.

FIG. 13 shows a process flow diagram to illustrate a process of howling detection according to one example. The process illustrated in FIG. 13 may, for example, take place in a howling detector 150 as illustrated in FIGS. 5, 6 and 8 .

In a first step ST1 a signal comprising audio data derived from a microphone is obtained and forms an input signal u(n). In a second step ST2 the log energy y(n) of the signal u(n) is computed from the amplitude a(n) using, e.g. equation (2). The signal y(n) is normalised at ST3 such that it is between 0 and 1, as given by equation 5. A line of best fit is derived at ST4 in order to fit a straight line to the normalised plot of the log energy. The resulting line or trend line is smoothed at ST5 using e.g. equation (9). Steps ST1 to ST5 inclusive can be considered to form part of a processing stage A which comprises the generation of a trend line.

In a decision-making stage B, a first decision D1 is made at ST6 in order to ascertain if the energy level y(n) is below a given threshold. If the energy level is below the threshold, a no-howling decision (or howling flag=0) is reached. If the energy level is determined to be at or above the energy threshold, a number of further decision making steps are carried out in order to determine a log-linearity metric. Specifically, at ST7 a goodness of fit measure of the computed trend line to the logarithm of the energy of the input signal over P samples is determined. The goodness of fit measure may be determined, for example, in accordance with equation (1) and the resulting linearity metric may be smoothed according to equation (8).

A second decision D2 is made at ST9. Specifically, the linearity metric is compared with a predetermined threshold wherein, if the linearity metric is below the predetermined threshold, indicating that a goodness of fit between the observed samples and the trend line returns high residual values, a no-howling decision is reached. However, if the linearity metric is at or above the predetermined threshold, a third decision-making process D3 is performed at ST10 in order to ascertain if the slope of the estimated trend line falls between the confidence intervals according to equation (11). If the slope does not fall between the confidence intervals, a no-howling decision is reached. However, if the slope does fall between the confidence intervals, a result of howling (or howling flag=1) is reached.

It will be appreciated that the indication of howling will be reached when the linearity metric exceeds a predetermined threshold and that the threshold may be set to detect circumstances of imminent howling (tending towards howling) as well as circumstances of actual howling. If the predetermined threshold is set to detect circumstances of imminent howling (in other words the linearity metric exceeds a lower threshold and/or may be tending towards a higher linearity metric threshold, thus indicating that the system becoming more unstable), the gain adjuster may be operable to adjust the gain by an amount proportional to the level of instability of the system, e.g. by a few decibels, rather than by a larger (typically fixed) amount that will interrupt/kill the noise control unit.

Thus, in a subsequent step a command may be issued to a gain adjustment unit in order to adjust the gain of an associated noise control circuit.

The possible actions taken in response to the process illustrated in FIG. 13 may be as follows:

No howling: no action taken.

Howling: Drop gain by fixed amount in order to interrupt ANC.

Imminent howling: reduce ANC gain slightly (e.g. 3-5 dB or by an amount derived from the slope of the fitted line (without any normalisation)).

It will be appreciated that a howling detection process according to one or more examples may be restricted to at least one defined frequency band—e.g. between f1 and f2. This may be achieved by filtering the input signal in order to eliminate certain frequencies or concentrate on certain frequency band(s). An advantage of this arraignment is that “howling-type” acoustic patterns can be filtered out. That is, audio sounds such as whistling or sirens which may otherwise give rise to a false positive indication of howling. Thus, an advantage of one or more of the present examples is that the occurrence of false positives is alleviated, resulting in a more robust system.

According to one or more examples, and as illustrated in FIG. 14 , howling detection may be carried out on the basis of a plurality of frequency bands. This may be referred to as multiband howling detection, wherein a plurality of filters 180 _(1 . . . M) allow the input signal u(n) to be tuned for a plurality M of frequency bands. Thus, the system may be considered to comprise a filterbank which carries out frequency band splitting of the input signal.

The howling detector 150 _(1 . . . M) is operable to determine a log-linearity metric for each frequency band. The slope and confidence measures may be readily combined in a combiner 190 by checking each band independently and then acting upon the channel with high confidence and slope closest to 1/(P−1). According to this arrangement, a decision of howling may be reached if howling is detected in any single frequency band.

The skilled person will recognise that some aspects of the above-described apparatus and methods may be embodied as processor control code, for example on a non-volatile carrier medium such as a disk, CD- or DVD-ROM, programmed memory such as read only memory (Firmware), or on a data carrier such as an optical or electrical signal carrier. For many applications examples of the invention will be implemented on a DSP (Digital Signal Processor), ASIC (Application Specific Integrated Circuit) or FPGA (Field Programmable Gate Array). Thus the code may comprise conventional program code or microcode or, for example code for setting up or controlling an ASIC or FPGA. The code may also comprise code for dynamically configuring re-configurable apparatus such as re-programmable logic gate arrays. Similarly the code may comprise code for a hardware description language such as Verilog TM or VHDL (Very high speed integrated circuit Hardware Description Language). As the skilled person will appreciate, the code may be distributed between a plurality of coupled components in communication with one another. Where appropriate, the examples may also be implemented using code running on a field-(re)programmable analogue array or similar device in order to configure analogue hardware.

Note that as used herein the term unit or module shall be used to refer to a functional unit or block which may be implemented at least partly by dedicated hardware components such as custom defined circuitry and/or at least partly be implemented by one or more software processors or appropriate code running on a suitable general purpose processor or the like. A unit may itself comprise other units, modules or functional units. A unit may be provided by multiple components or sub-units which need not be co-located and could be provided on different integrated circuits and/or running on different processors.

Examples may be implemented in a host device, especially a portable and/or battery powered host device such as a mobile computing device for example a laptop or tablet computer, a games console, a remote control device, a home automation controller or a domestic appliance including a smart home device a domestic temperature or lighting control system, a toy, a machine such as a robot, an audio player, a video player, or a mobile telephone for example a smartphone.

Examples of the disclosure may be provided according to any one of the following numbered statements:

1. A processing module for a noise control circuit, the processing module comprising:

-   -   a howling detector configured to receive an input signal and to         determine a linearity metric based on the input signal, the         linearity metric comprising a measure of the linearity of a         logarithmic representation of the energy of the input signal;         and     -   a gain adjuster configured to adjust the gain of a noise control         unit.

2. The processing module of statement 1, wherein the linearity metric is derived by computing a difference between the logarithm of the energy of the signal and a trend line, the trend line being a straight line representation of the logarithm of the energy of the signal.

3. The processing module of statement 2, wherein coefficients of the trend line are derived by a fitting a line to the logarithmic representation of the energy of the input signal.

4. The processing module of statement 3, wherein fitting a line to the logarithmic representation of the energy of the input signal comprises performing a least squares computation which minimises a difference between the trend line and the logarithmic representation of the energy of the input signal.

5. The processing module of any one of statements 2 to 4, wherein the linearity metric is determined according to a goodness of fit measure of the computed trend line to the logarithm of the energy of the input signal over P samples.

6. The processing module of statement 5, wherein the goodness of fit measure r(n) is represented by:

${r(n)} = \sqrt{\sum\limits_{i = 0}^{P - 1}\left( {{y\left( {n - i} \right)} - \left( {{{\hat{\beta}}_{0}(n)} + {\left( {P - i} \right) \cdot {{\hat{\beta}}_{1}(n)}}} \right)} \right)^{2}}$ where y(n) is the log energy of the input signal, {circumflex over (β)}₀(n) is the bias of the trend line and {circumflex over (β)}₁(n) is the slope of the trend line.

7. The processing module of any preceding statement, wherein the howling detection unit is further configured to issue a command to the gain adjuster in order to reduce the gain if the linearity metric exceeds a predetermined threshold.

8. The processing module of any preceding statement, wherein the gain adjuster is configured to adjust the gain by a fixed amount.

9. The processing module of any preceding statement, wherein the howling detection unit is further configured to estimate a maximum stable gain of the noise control unit.

10. The processing module of statement 9 when appended directly or indirectly to statement 2, wherein the maximum stable gain is proportional to the slope of the trend line.

11. The processing module of any preceding statement, wherein the howling detection unit is configured to determine, based on the determined linearity metric, if howling is likely or imminent.

12. An audio processing system comprising a processing module of any preceding statement and further comprising a noise control unit for generating a noise control signal based on a reference input signal which represents a sound detected by a reference microphone.

13. The audio processing system of statement 12, wherein the processing module is connected to a speaker and wherein the speaker generates an anti-noise signal based on the noise control signal in order to cancel or at least reduce the noise detected by a reference microphone.

14. The audio processing system of statement 13, wherein the howling detector is provided in parallel with the noise control unit.

15. The audio processing system of any of statements 12-14, further comprising a filter configured to filter out one or more frequencies or frequency bands of the input signal.

16. The audio processing system of any of statements 12-15, further comprising a filterbank configured to split the input signal into a plurality of frequency bands, wherein the howling detector is operable to determine a linearity metric for each frequency band.

17. A processing circuit for a noise control module comprising:

a gain adjustment mechanism configured to adjust the gain of a noise control circuit if a plot of the energy of an input signal in the log domain becomes linear or tends towards linearity.

18. The processing circuit of statement 17, wherein the gain of the noise control circuit is adjusted by an amount that relates to a difference between a measure of the linearity of the plot of the energy of the input signal in the log domain and a target linearity metric.

19. The processing module ofany one of statements 1 to 15, in the form of a monolithic integrated circuit.

20. A device comprising a processing module according to any one of statements 1 to 15, or a processing circuit according to statement 17 or 18, wherein the device comprises a mobile telephone, headphone, acoustic noise cancelling headphones, a smart watch, an audio player, a video player, a mobile computing platform, a games device, a remote controller device, a toy, a machine, or a home automation controller, a domestic appliance or other portable device.

21. A method of processing an audio signal comprising:

-   -   determining a linearity metric based on an input signal, the         linearity metric comprising a measure of the linearity of a         logarithmic representation of the energy of the input signal;         and     -   adjusting the gain of the noise control unit if the linearity         metric exceeds a predetermined threshold.

22. The method of processing an audio signal of statement 21, wherein the step of determining a linearity metric comprises computing a difference between the logarithm of the energy of the signal and a trend line, the trend line being a straight line representation of the logarithm of the energy of the signal.

It should be noted that the above-mentioned examples illustrate rather than limit the invention, and that those skilled in the art will be able to design many alternative examples without departing from the scope of the appended claims. The word “comprising” does not exclude the presence of elements or steps other than those listed in a claim, “a” or “an” does not exclude a plurality, and a single feature or other unit may fulfil the functions of several units recited in the claims. Any reference numerals or labels in the claims shall not be construed so as to limit their scope. 

The invention claimed is:
 1. A howling detector configured to: receive an input audio signal; determine a parameter of a portion of the input audio signal; determine a measure of linearity of the parameter; and compare the measure of linearity to a predetermined threshold to determine whether howling is likely or imminent on the basis of a parameter of a portion of an input signal.
 2. The howling detector of claim 1 wherein the howling detector is configured to measure the parameter of the portion of the input signal.
 3. The howling detector of claim 1 wherein the howling detector is configured to determine whether howling is likely or imminent on the basis of the linearity of the parameter of the portion of the input signal.
 4. The howling detector of claim 3 wherein the howling detector is configured to determine the linearity of the parameter of the portion of the input signal.
 5. The howling detector of claim 1 wherein the parameter is proportional to the energy of the portion of the input signal.
 6. The howling detector of claim 5 wherein the parameter is proportional to the logarithm of the energy of the portion of the input signal.
 7. The howling detector of claim 6 wherein the parameter of the portion of the input signal is proportional to the difference between the logarithm of the energy of the portion of the input signal and a trend line, wherein the trend line comprises a straight line representation of the logarithm of the energy of the portion of the input signal.
 8. The howling detector of claim 7 wherein the parameter of the portion of the input signal is proportional to the residuals of the trend line.
 9. The howling detector of claim 7 wherein the parameter of the portion of the input signal is proportional to a goodness of fit measure of the trend line over P samples.
 10. The howling detector of claim 9 wherein the goodness of fit measure, r(n), is represented by: ${{r(n)} = \sqrt{\sum\limits_{i = 0}^{P - 1}\left( {{y\left( {n - i} \right)} - \left( {{{\hat{\beta}}_{0}(n)} + {\left( {P - i} \right) \cdot {{\hat{\beta}}_{1}(n)}}} \right)} \right)^{2}}},$ where y(n) is the log energy of the input signal, {circumflex over (β)}₀(n) is the bias of the trend line and {circumflex over (β)}₁(n) is the slope of the trend line.
 11. The howling detector of claim 7 wherein the parameter is proportional to the slope of the trend line.
 12. The howling detector of claim 5 wherein the parameter is proportional to a derivative of the energy of the portion of the input signal.
 13. The howling detector of claim 6 wherein the parameter is proportional to a derivative of the logarithm of the energy of the portion of the input signal.
 14. A gain adjuster configured to adjust the gain of a noise control circuit based on a command issued by a howling detector indicating that a measure of linearity of a parameter of a portion of an input audio signal exceeds or falls below a predetermined threshold which indicates whether howling is likely or imminent.
 15. The gain adjuster of claim 14 wherein the gain adjuster is configured to adjust the gain of a noise control circuit based on the linearity of the parameter of the portion of the input signal.
 16. The gain adjuster of claim 14 wherein the parameter is proportional to the energy of the portion of the input signal.
 17. The gain adjuster of claim 14 wherein the parameter is proportional to the logarithm of the energy of the portion of the input signal. 